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The  purpose  of  this  research  was  to  determine  if  favorable  interference 
effects  could  be  obtained  from  a wing-rotor  interaction  system.  It  was 
theorized  by  the  fifth  author,  J.  J.  Murray,  that  when  this  system  is  in 
forward  flight  the  oncoming  stream  would  sweep  the  wake  from  the  rotor  over 
the  top  of  the  wing  which  would  increase  the  velocity  on  the  wing.  The 
additional  velocity  passing  over  the  top  of  the  wing  would  decrease  the 
pressure  and  increase  the  lift.  This  additional  velocity  would  also  delay 
boundary  layer  separation,  increase  the  stall  angle  of  attack,  and  decrease 
the  drag. 

This  final  report  covers  the  research  performed  from  September  1,  1973 
to  December  31,  1977.  The  first  two  years  was  devoted  to  an  experimental 
investigation  in  the  North  Carolina  State  University  Subsonic  Wind  Tunnel. 

A wing  which  spanned  half  of  the  test-section  width  was  mounted  in  a canti- 
lever fashion  from  a side  wall,  and  the  proprotor  was  positioned  above  the 
wing  tip  on  a shaft  which  ran  through  the  ceiling.  The  wing  and  proprotor 
were  instrumented  separately  to  measure  their  lift  and  drag,  and  one  wing 
model  was  constructed  to  measure  the  pressure  distribution.  The  test  re- 
sults showed  significant  increases  in  lift  and  decrease  in  drag  for  both 
the  wing  and  the  proprotor  due  to  favorable  interference  effects.  Measured 
pressure  distributions  on  the  wing  revealed  that  the  inpingement  of  the 
proprotor 's  wake  delayed  separation  and  thus  infcreased  the  stall  angle  of 
the  wing.  A complete  description  of  this  investigation  is  given  in  the 
Master's  thesis  (ref.  1)  of  S.  L.  Griffith,  fourth  author.  A condensed 
version  was  presented  at  the  AIAA  3rd  Atmospheric  Flight  Mechanics  Conference 
held  in  Arlington,  Texas,  June  7-9,  1976  and  appears  in  the  proceedings  of 
that  conference  (ref.  2). 
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The  experimental  Investigations  showed  that  the  concept  of  favorable 
interference  with  the  wing-rotor  system  was  feasible  and  prompted  an  analyti- 
cal investigation.  In  order  to  calculate  the  combined  flow  field  of  the 
rotor  and  the  wing,  it  is  necessary  to  have  a method  which  will  calculate 
subsonic  wing  properties  accurately  and  expediously.  The  last  two  years 
of  this  research  was  devoted  to  this  task. 

The  first  author  (ref.  3)  developed  a new  lifting  surface  theory  for 
rectangular  wings.  The  solution  can  also  be  interpreted  as  a vortex  lattice 
method;  however,  the  optimum  locations  of  horseshoe  vortices  and  control 
points  are  determined  in  the  solution.  It  properly  locates  the  tip  vortices 
to  be  inset  from  the  wing  tips  which  other  investigators  found  empirically. 

The  method  of  solution  properly  accounts  for  the  leading-edge  singularity, 

Cauchy  singularity,  and  Kutta  condition.  Results  were  found  to  be  more 
accurate  and  converge  faster  than  conventional  vortex  lattice  methods.  In 
addition,  the  accuracy  was  found  to  be  comparable  to  other  lifting  surface 
theories,  but  with  much  smaller  computational  times.  This  method  was  pre- 
sented at  the  Vortex  Lattice  Utilization  Workshop  held  at  NASA  Langley 
Research  Center,  May  17-18,  1976.  A description  of  the  technique  used  to 
extract  the  leading-edge  suction  parameter  from  the  solution  is  given  in 
reference  4,  and  it  was  also  presented  at  13th  Annual  Meeting  of  the  Society 
of  Engineering  Science,  Hampton,  Virginia,  November  1-3,  1976. 

The  first  three  authors  extended  the  lifting  surface  theory  to  skewed 
and  swept  wings  (ref.  5 attached  to  this  report).  The  extension  to  skewed 
wings  was  straightforward;  however,  the  midspan  point  of  swept  wings  presented 
a problem.  A "kink"  in  the  spanwise  vortex  occurs  there,  and  the  theory  must 
be  modified  in  a fashion  similar  to  that  used  by  Multhopp  in  order  to  get  a 
good  solution.  The  spanwise  variation  of  vorticity  was  calculated  by  two 
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methods,  a Fourier  sine  series  and  a direct  solution  for  the  vorticity  at 
specific  spanwise  locations.  Both  methods  yielded  essentially  the  same 
results  for  rectangular  wings,  but  the  Fourier  sine  series  gave  oscillatory 
results  for  the  swept  wing.  In  general,  the  direct  method  gave  satisfactory 
results  which  compared  well  with  other  methods.  A condensed  version  of 
reference  5 will  be  submitted  to  the  Journal  of  Aircraft. 

RECOMMENDATIONS 

It  is  recommended  that  the  analytical  investigation  be  extended  to  in- 
clude a rotor  passing  over  the  wing.  A more  accurate  mathematical  model  of 
the  wake  should  also  be  included. 
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1 . SUMMARY 


A new  method  1s  developed  for  solving  the  lifting  surface  equation 
for  thin  wings.  The  downwash  equation  is  transformed  into  double  inte- 
grals Involving  Cauchy-type  singularities  in  the  chordwise  and  spanwise 
directions.  A technique  developed  by  Lan  for  airfoil  theory  is  used  to 
reduce  both  Integrals  to  a double  summation.  This  method  properly  ac- 
counts for  the  leading-edge  singularity,  Cauchy  singularity,  and  Kutta 
condition.  The  circulation  is  determined  by  solving  the  matrix  equation 
formed  by  requiring  the  downwash  equatiori  to  satisfy  the  tangent-flow 
boundary  condition  at  specific  control  points  on  the  wing.  Two  methods 
were  used  for  the  spanwise  variation  of  the  circulation:  a Fourier-sine 
series  and  a direct  matrix  solution  which  determines  the  circulation  at 
specific  spanwise  locations.  These  two  methods  were  found  to  yield  essen- 
tially the  same  results  for  rectangular  and  skewed  wings.  However,  the 
direct  matrix  method  gave  more  satisfactory  results  for  swept  wings.  The 
differences  are  attributed  to  the  infinite  downwash  caused  by  the  "kink" 
in  the  spanwise  vortices  at  the  midspan  of  swept  wings.  The  solutions 
generally  compared  well  with  other  lifting-surface  theories,  but  with 
much  smaller  computational  times,  and  the  method  was  found  to  be  more 
accurate  and  converge  faster  than  conventional  vortex  lattice  methods. 


The  motivation  for  the  research  described  in  this  report  stemmed  from 
an  experimental  investigation  of  a wing-rotor  interaction  system  for  heli- 
copters. A complete  description  of  this  investigation  is  given  in  reference  1, 
and  a condensed  version  appears  in  reference  2.  Subsonic  wind  tunnel  tests 
were  performed  to  determine  the  interference  effects  of  a wing  and  proprotor 
on  their  lift  and  drag  characteristics.  A wing  which  spanned  half  of  the 
test-section  width  was  mounted  in  a cantilever  fashion  from  a side  wall,  and 
the  proprotor  was  positioned  above  the  wing  tip  on  a shaft  which  ran  through 
the  ceiling.  The  wing  and  proprotor  were  instrumented  separately  to  measure 
their  lift  and  drag,  and  one  wing  model  was  constructed  to  measure  the 
pressure  distribution.  The  test  results  showed  significant  increases  in  lift 
and  decrease  in  drag  for  both  the  wing  and  the  proprotor  due  to  favorable 
interference  effects.  Measured  pressure  distributions  on  the  wing  revealed 
that  the  impingement  of  the  proprotor's  wake  delayed  separation  and  thus  in- 
creased the  stall  angle  of  the  wing.  These  features  made  the  concept  attrac- 
tive for  future  design  considerations  and  prompted  the  present  analytical 
investigation  of  subsonic  wings. 

An  analytical  investigation  of  a complete  wing  rotor  flow  field  requires 
computational  methods  for  both  the  wing  and  the  rotor.  In  order  to  calculate 
the  combined  flow  field  efficiently,  it  is  necessary  to  have  a method  which 
will  calculate  subsonic  wing  properties  accurately  and  expediously.  It  is 
the  purpose  of  this  report  to  develop  such  a method. 

Aerodynamic  characteristics  of  wings  at  subsonic  speeds  are  generally 
calculated  by  vortex  lattice  or  lifting  surface  methods.  In  the  conventional 
vortex  lattice  method  (VLM)  the  planform  of  the  wing  is  divided  into  a number 
of  elemental  panels,  and  a horseshoe  vortex  is  placed  at  the  local  quarter- 
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chord  of  each  panel.  The  boundary  condition  Is  satisfied  at  the  local  three- 
quarter  chord  of  each  panel  by  requiring  the  flow  to  be  tangent  to  the  surface 
there.  The  strengths  of  the  horseshoe  vortices  are  determined  by  solving  the 
matrix  equation  formed  from  the  boundary-condition  equations.  Then  the  aero- 
dynamic characteristics,  such  as  lift,  Induced  drag,  pitching  moment,  and 
rolling  moment  are  calculated  by  summing  the  results  from  each  panel.  A 
complete  description  of  the  vortex  lattice  method  Is  given  In  reference  3. 
Vortex- lattice  methods  are  simpller  and  easier  to  apply  to  complex  configura- 
tions than  lifting  surface  methods,  but  they  are  generally  less  accurate 
(ref.  4). 

In  lifting  surface  theories,  a continuous  loading  or  circulation  In 
both  the  spanwise  and  chordwise  directions  Is  used  to  determine  the  aero- 
dynamic characteristics  of  the  wing.  Here,  singular  Integral  equations  must 
be  solved  which  cause  numerical  difficulties.  The  load  distribution  Is 
generally  represented  by  a double  series  and  the  coefficients  of  these 
series  are  determined  by  solving  the  matrix  equation  which  results  from 
applying  the  tangent-flow  boundary  condition  at  specific  points,  called  con- 
trol points.  Lifting-surface  theories  are  discussed  In  reference  5. 

Lan  (ref.  4)  developed  an  Ingenious  method  for  thin,  two-dimensional 
airfoils  by  using  the  midpoint  trapezoidal  rule  and  the  theory  of  Chebychev 
polynomials  to  reduce  the  downwash  Integral  to  a finite  sum.  This  method  Is 
simple  to  apply  and  gives  an  accurate  solution  to  the  airfoil  downwash  equa- 
tion. A more  detailed  description  Is  given  below.  Lan  also  developed  a 
quasi-vortex  lattice  method  for  finite  wings  by  using  his  airfoil  method 
for  a continuous  chordwise  vortex  distribution  In  the  spanwise  direction. 

The  results  showed  an  Improvement  over  the  conventional  vortex  lattice 
method.  An  additional  Improvement  was  found  empirically  by  Hough  (ref.  6) 
when  the  lattice  was  Inset  1/4  of  a lattice  width  at  the  wing  tips. 
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DeJarnette  (ref.  7)  developed  a new  lifting  surface  theory  for  rectan- 
gular wings  which  uses  Lan's  continuous  chordwise  vortex  distribution  but, 
unlike  Lan,  a continuous  spanwise  vortex  distribution  also.  Although  the 
vortex  distributions  are  continuous,  the  method  Is  easily  Interpreted  as  a 
vortex  lattice  method  In  which  the  arrangement  of  horseshoe  vortices  and 
control  points  are  determined  from  the  solution  rather  than  chosen  like  the 
conventional  vortex  lattice  method.  It  properly  locates  the  tip  vortices  to 
be  Inset  from  the  wing  tips  like  Hough  (ref.  6)  found  empirically.  The  re- 
sults were  found  to  compare  very  closely  with  other  lifting  surface  theories, 
but  the  computational  time  was  only  about  1/10  of  the  other  methods.  The 
results  were  more  accurate  and  converged  much  more  rapidly  than  the  vortex 
lattice  method  of  reference  3.  A description  of  the  method  used  for  ex- 
tracting the  leading-edge  suction  parameter  from  the  solution  Is  given  In 
reference  8. 

The  lifting  surface  theory  given  In  references  7 and  8 apply  the  down- 
wash  Integral  In  the  form  Involving  Cauchy  singularities  rather  than  the 
Mangier  singularity  form  which  Is  usually  used  In  lifting  surface  theories. 
Consequently,  difficulties  are  encountered  when  the  same  method  Is  applied  to 
non-rectangular  planforms.  This  report  develops  a generalized  lifting  sur- 
face theory  which  can  be  applied  to  wings  which  are  skewed,  swept,  or  tapered, 
yet  It  retains  the  same  accuracy  and  simplicity  as  the  method  given  In 
references  7 and  8. 


3.  SYMBOLS 


parameter  defined 

by 

eq. 

(7.13) 

parameter  defined 

by 

eq. 

(7.26) 

unknowns  In  the  downwash 

matrix 

parameter  defined 

by 

eq. 

(6.10) 

parameter  defined 

by 

eq. 

(6.14) 

parameter  defined 

by 

eq. 

(8.12) 

parameter  defined 

by 

eq. 

(B-3) 

parameter  defined 

by 

eq. 

(8.19) 

aspect  ratio 

wing  span 

parameter  defined 

by 

eq. 

(7.14) 

parameter  defined 

by 

eq. 

(7.27) 

parameter  defined 

by 

eq. 

(6.11) 

parameter  defined 

by 

eq. 

(6.15) 

parameter  defined 

by 

eq. 

(8.13) 

parameter  defined 

by 

eq. 

(B-4) 

parameter  defined 

by 

eq. 

(8.20) 

wing  chord 
mean  chord  of  wing 
wing  root  chord 
wing  tip  chord 
left  wing  tip  chord 
right  wing  tip  chord 
sectional  lift  coefficient 
wing  lift  coefficient 
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Cl 

(X 

lift  curve  slope,  per  radian 

Cm 

sectional  moment  coefficient 

Cm 

wing  pitching  moment  coefficient,  about  leading  edge 
at  midspan 

Cs 

leading  edge  suction  parameter 

d 

parameter  defined  by  eq.  (7.10) 

d* 

parameter  defined  by  eq.  (7.23) 

e 

parameter  defined  by  eq.  (7.11) 

e' 

parameter  defined  by  eq.  (7.24) 

f 

parameter  defined  by  ew.  (7.12) 

f 

parameter  defined  by  eq.  (7.25) 

G 

parameter  defined  by  eq.  (5.8) 

parameter  defined  by  eq.  (6.13) 

‘•ijkJi 

parameter  defined  by  eq.  (8.18) 

\c 

^Sjka 

parameter  defined  by  eq.  (8.16) 

parameter  defined  by  eq.  (8.11) 

parameter  defined  by  eq.  (B-2) 

^RC 

parameter  defined  by  eq.  (8.9) 

•"le 

tan 

"*TE 

tan 

Am 

tan  A^^  - tan  Aj^^ 

M 

number  of  spanwise  vortices 

n 

summational  Integer 

N 

number  of  chordwise  vortices  and  control 

points 

S 

wing  planform  area 

u,v 

parameters  for  integration  by  parts,  see 

eq.  (7.4) 

VLM 

vortex  lattice  method 

V 


freestream  velocity 

doMTiMash  velocity,  positive  upward 

downwash  velocity  Induced  by  the  vortices  on  left  wing  panel 
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I 

1 


X 


y 

Y 

Y' 


a 


Y 

r 


n 


0 


X 

A 

LE 

A 

TE 

♦ 


downwash  velocity  Induced  by  the  left  vortex  filament  at  the 
wing  midspan 

downwash  velocity  Induced  by  the  vortices  on  right  wing  panel 

downwash  velocity  Induced  by  the  right  vortex  filament  at 
the  wing  midspan 

chordwise  coordinate  measured  from  leading  edge  In  direction 
of  V 

00 

parameter  defined  by  eq.  (6.9) 
wing  aerodynamic  center  location 

spanwise  coordinate  on  right  wing  panel,  positive  to  the  right 

spanwise  coordinate  on  left  wing  panel 

parameter  defined  by  eq.  (7.8) 

parameter  defined  by  eq.  (7.21) 

parameter  defined  by  eq.  (7.9) 

parameter  defined  by  eq.  (7.22) 

vertical  coordinate  of  mean  camber  line 

angle  of  attack 

non-dimensional  circulation  per  unit  chord 
circulation 

parameter  defined  by  eq.  (9.18) 

non-dimensional  chordwise  coordinate,  see  eq.  (6.2) 

transformed  chordwise  coordinate 

wing  taper  ratio 

leading  edge  sweep 

trailing  edge  sweep 

transformed  spanwise  coordinate  on  right  wing  panel 
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‘ i '■ 


Subscripts 


transformed  spanwise  coordinate  on  left  wing  panel 

transformed  spanwise  coordinate  for  entire  span,  see  eqs.  (8.2) 
and  (8.3) 


chordwise  control  point,  see  eq.  (4.7) 
spanwise  control  point,  see  eq.  (5.14) 
chordwise  vortex  point,  see  eq.  (4.6) 


spanwise  vortex  point,  see  eq.  (5.13) 
leading  edge  of  wing 
trailing  edge  of  wing 
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4.  LAN'S  TWO-DIMENSIONAL  THEORY 


For  thin  airfoils,  the  downwash  equation  Is 


1 } ‘*'‘1 
^ - it  J — ^ 


(4.1) 


The  Integral  on  the  right  side  Is  of  the  Cauchy  type,  and  the  Integrand 
has  a singularity  at  x-j  “ Xj  • Transform  the  x coordinate  by 


x/c  * (1  - cos  e)/2 

and  use  the  following  result  from  airfoil  theory  (ref.  9) 


(4.2) 


f 

J cos  0^  - 


cos  e 


(4.3) 


to  write  eq.  (4.1)  as 


^ J Y(e|)  sin  0^  d0^  j J [Y(e^)  sin  0^  - y(0)  sin  0]  d0^ 

'^l  ~ ” lir  J cos  0,  - cos  0 “ ” lir  J COS  0,  - COS  0 

0 0 

Now  the  singularity  In  the  Integrand  has  been  removed  because  the  numerator 
approaches  zero  at  the  same  position  (0^  = 0)  that  the  denominator 
approaches  zero. 

Lan  (ref.  4)  used  the  theory  of  Chebychev  polynomials  to  show  that 


cos  0|^  - cos  0^ 


when  the  vortex  positions  are 


[2k  - l)ir 


-N^  for  1=0 
0 for  1 / 0,  N 
N^  for  1 = N 


, k = 1 , . . . , N 


(4.5) 


(4.6) 
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and  the  control  points  are  located  at 


Then  the  Integral  In  eq.  (4.4)  can  be  reduced  to  a finite  sum  by  using  the 
midpoint  trapezoidal  rule  (ref.  4)  and  eq.  (4.5)  to  obtain 


1 . ? T|,  S)|.  - T,  Sin  e, 

Zn  N COS  0|^  - COS  0^ 


1 


Yk 


sin  0, 


cos  0|^  - cos  0^ 


U 11m  y(0)  sin  0 
^ 0^ 


7 11m  y(e)  sin  0 
^ 0-Hr 


However, 


1 = 0 

1 / 0,N  (4.8) 
1 = N 


11m  y(e)  sin  0 = 4 Cf.  (4.9) 

0->O  ^ 


where  Is  the  leading-edge  suction  parameter  and  since  the  Kutta  condition 
requires  that  ^(Tr)  = 0, 


11m  y(0)  sin  0=0.  (4.10) 

0-»ir 

Unlike  the  conventional  vortex  lattice  method,  the  Cauchy  singularity, 
leading-edge  square-root  singularity  and  the  Kutta  condition  are  properly 
accounted  for  In  this  method.  Equation  (4.8)  can  be  solved  with  1 f 0 to 
obtain  the  N values  of  Y|^»  and  then  the  leading-edge  suction  parameter  can 
be  computed  by  using  eq.  (4.8)  with  1 = 0 (control  point  at  the  leading 
edge).  Figure  1 Illustrates  the  positions  of  the  vortices  and  control 
points  by  the  "semi -circle  method"  for  N * 2.  With  only  one  vortex 
(N  • 1)  on  a flat  plate,  the  exact  lift  and  leading-edge  suction  are  ob- 
tained, and  the  Kutta  condition  Is  satisfied.  With  two  or  more  vortices 
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the  exact  pitching  moment  Is  obtained  In  addition  to  the  above  properties 
and  the  calculated  Y|^  are  exact.  It  can  be  shown  that  the  remarkable 
accuracy  of  this  method  Is  due  to  eq.  (4.5),  which  is  similar  to  the  Inte- 
gral result  given  by  eq.  (4.3)  and  used  In  exact  thin  airfoil  theory. 

5.  PRANDTL'S  LIFTING-LINE  THEORY 

Before  developing  the  lifting  surface  equations  for  skewed  and  swept 
wings,  the  technique  for  performing  the  spanwise  Integration  will  be 
tested  with  Prandtl's  lifting-line  equation.  The  lifting-line  equation 
Is  given  by  ref.  9 as 


r = ttV„c 


a 


dy^  (y-y^ ) 


(5.1) 


This  equation  has  a Cauchy  type  integral  on  the  right  side,  and  therefore 
the  method  used  by  Lan  for  the  airfoil  is  applicable.  In  reference  7 the 
spanwise  coordinate  was  transformed  by  the  "semi-circle"  method 


y = - I cos  <(,  (5.2) 

and  the  circulation  r((|>)  was  represented  by  Multhopp's  interpolation  formula 

which  is  a Fourier  sine  series.  Then  Lan's  airfoil  technique  was  used  to 

solve  equation  (5.1).  The  results  were  found  to  be  identical  to  the  classi- 
cal solution  of  Prandtl's  lifting-line  equation  (ref.  9). 

When  using  the  Fourier  sine  series  for  r in  eq.  (5.1),  the  series 
must  be  differentiated  and  the  accuracy  of  the  resulting  series  is  ques- 
tionable. Instead  of  using  a series  for  r,  values  of  the  derivative  of  r 
at  discrete  spanwise  locations  will  be  calculated  here.  When  treating 
swept  wings  In  Chapter  7,  It  will  be  necessary  to  divide  the  spanwise 
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I 1 

K ^ 


Integration  range  Into  a left-wing  panel  and  a right-wing  panel.  There- 
fore, for  the  analysis  here  the  spanwise  coordinate  will  be  transformed 
by  applying  the  "semi -circle"  method  to  each  wing  panel  separately.  As 
shown  In  Appendix  A,  the  transformation  given  by  eq.  (5.2)  Is  replaced 


by  the  following: 


y = j (1  - cos  (|»)  for  y > 0 (right-wing  panel) 


y = - ^ (1  - cos  i)  for  y < 0 (left-wing  panel) 


(5.3) 


(5.4) 


The  downwash  Integral  In  eq.  (5.1)  can  be  divided  Into  a contribution 
from  the  left-wing  panel  and  right-wing  panel. 

w.  * Wp  + w, 

and  for  yj  > 0, 


b/2  j TT  j. 

w = ^ f dr  , - 1 4 f dr  ^^1 

J (y  .-y, ) 4TrV  b J d4  (cos  4,,  - cos 
j n ' J ' “ n ' J 


(5.6) 


0 0 

- f dr  ^^1  . 1 4 f dr 

■ 4irv„  J ^ (yj-y^)  ■ ^ H J 


COS  ((ij  - cos 


(5.7) 


Note  that  eq.  (5.6)  has  the  Cauchy-type  singularity  but  eq.  (5.7)  does  not 
have  a singularity  In  the  Integrand  because  y^  > 0 and  the  Integral  covers 
the  left-wing  panel  where  y-j  < 0.  With 


(5.8) 


a wing  symmetrically  loaded  about  the  midspan  line  (y  » 0)  will  have 


G(^  * GCifr) 

and  thus  equation  (5.7)  can  be  written  as 


(5.9) 
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IT 


G('>l)  d(>^ 

- cos  ♦j  - cos 


(5.10) 


Now  replace  eqs.  (5.6)  and  (5.10)  with  the  midpoint  trapezoidal  rule 
summation  to  get 


•'r,  * ' FBR 


1 w 

TR  I. 


[cos  4.^  - cos  (|)j, 


(5.11) 


"Lj  ■ *3*  J, 


cos  - cos 


(5.12) 


These  two  equations  represent  the  downwash  due  to  trailing  vortices  located 


7 - * - (2t  - Dir  -1 

*i,  ~ ZM  ’ t - 1,  ....  M 


(5.13) 


with  control  points  located  at 


^ “ 0*  •••> 


(5.14) 


Since  the  downwash  Is  expressed  In  terms  of  the  unknown  = (dr/d4>)j^. 
It  Is  necessary  to  express  r In  terms  of  G before  a solution  to  eq.  (5.1) 
can  be  obtained.  Noting  that  r = 0 at  the  wing  tips,  replace  r by 


f 

J ^ 


(5.15) 
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and  approximate  the  Integral  by  the  midpoint  trapezoidal  rule  to  get 


(5.16) 


Finally,  substitute  eqs.  (5.11),  (5.12),  and  (5.16)  Into  (5.1)  to  obtain 
. M M 


r I I 

t=3+i  ^ t=i  ^ 


1 


1 


2 - cos  ♦ - cos  (ji  . " cos  (|>  - cos  (}>  . 

X.  J J 


-oV  Mb 


(5.17) 


Apply  eq.  (5.17)  at  the  spanwise  locations  <)>.(j  = 0,1 M-1)  to  obtain 

a matrix  equation  for  the  M values  of  Gj^(i!-  = 1,  ...,  M).  Then  the  aero- 
dynamic coefficients  can  be  calculated  by  reducing  the  spanwise  Integrals 
to  a finite  sum  through  the  midpoint  trapezoidal  rule.  For  example,  the 
lift  coefficient  Is  given  by 


C 


L 


2V  V b 

00  J , 00 

~ j ^ *^1  " “5“ 

-b/2 


TT 

I T sin  ij>i  d(j)^ 
0 


(5.18) 


Integrate  the  right  side  by  parts  to  get 

Vcob 

= “T 

Replace  the  right  side  by  the  midpoint  trapezoidal  rule  to  get 


w 

^.^=0  ■ I ^ ^0^  *^1  ‘*'•’1 
0 


(5.19) 


M 

I 

)l=l 


r 

■ M 


M 


(5.20) 


Other  aerodynamic  coefficients  can  be  calculated  In  a similar  fashion. 


r r 


6.  LIFTING  SURFACE  THEORY  FOR  RECTANGULAR  AND  SKEWED  WINGS 


The  downwash  equation  from  lifting  surface  theory  (ref.  10),  In 
the  Mangier  Integral  form.  Is 


w(x.y)  = ^ 


S /(x-x^)^+{y-y^)' 


dx^  dy^ 


(6.1) 


Eq.  (6.1)  can  be  generalized  to  Include  skewed  wings  with  taper  by  trans- 
forming the  Independent  variables  from  x,y  to  n,y  where 


i(x,y) 


x-x,c(y) 


(6.2) 


Is  the  nondimenslonal  chordwise  variable  (0  < n < 1).  Thus, 


X = ^ F(y)  + n c(y) 


(6.3) 


and  eq.  (6.1)  becomes 


, . 1 rfY(ni.yi)  ^ x-x. 


/(x-x^)^+(y-y^)^ 


where  x = x(n,y) 


(6.4) 


Equation  (6.4)  Is  not  a Cauchy-type  Integral  due  to  the  term  (y-y^)^ 
In  the  denominator.  In  order  to  apply  Lan's  two-dimensional  method  for 
the  chordwise  Integration  and  the  lifting-line  method  for  the  spanwise 
Integration,  eq.  (6.4)  must  be  transformed  Into  a Cauchy-type  Integral 
form.  To  make  this  transformation.  It  Is  necessary  to  define  the  geometry 
and  coordinate  system  of  the  skewed  wing  as  shown  In  figure  2,  and  then 
Integrate  eq.  (6.4)  by  parts  In  the  spanwise  direction. 
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Defining  m|^£  = tan  Aj^^,  niy^  = tan  Ay^,  and  Am  = my^  - mj^^,  the  geometry 
Is  described  by 

= y "Jle  . ^TE  " y "h-E  (6.5) 

C(y)  = Xy£  - Xle 

= yAm  + Cp  (6.6) 

where  Cp  = ^ , Am  = j , and  X = . 

Substitute  eqs.  (6.5)  and  (6.6)  into  eq.  (6.3)  to  get 

X = y miE  + n[yAm  + Cj^]  (6.7) 

Substitute  eq.  (6.7)  into  eq.  (6.4)  and  integrate  by  parts  in  the  spanwise 
direction.  The  downwash  equation  then  becomes 


b/2 

l9yi  J 

r- 

w(n.y)  ■ - ^ J 

J 

1 

y-y^ 

^1 

( 

-b/2 

where 


X=  (A^  + 1)  (y-y^)^  + 2Ag  (y-y^)  + 


(6.8) 


(6.9) 


As  = + n-iAm  (6.10) 

Bg  = (n-ni)(C|^  + yAm)  (6.11) 

Equation  (6.8)  contains  Cauchy-type  integrals,  and  it  is  therefore 
in  the  form  to  apply  Lan's  technique.  With  y = -b/2  cos  (|>  and  n = (1-cos  e)/2, 
replace  both  integrals  in  eq.  (6.8)  with  the  midpoint  trapezoidal -rule  sum 
to  get 
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’ I. 


I i. 


M N 


^ -1  nn  f?  ? 

'*1  ,j  4irAR  M H 3*^  COS  (|»j-COS  <j) . 

' k j®"  " 


(6.12) 


where 


^ ^ /(Agi^+DCcos  »,^-cos  + 2Ag|^  Bgl^  (cos  - cos  »j)  ^^ 


''sk  ■ '"le  * f <’  - *k> 

Bsk  = (cos  - cos  0^)(^  - Y '*’j^ 


(6.13) 


(6.14) 


(6.15) 


The  control  points  are  located  at 


9^*-^,i  = l,...,N 


M 


, j = 1 M-1 


(chordwise) 


(s panwise) 


(6.16) 


(6.17) 


and  the  integration  or  vortex  points  are  located  at 


®k  “ •••»  N (chordwise) 

, £ = 1 M (spanwise) 


(6.18) 


(6.19) 


The  spanwise  variation  of  Y|^(<<>)  at  the  chordwise  point  0|^  is  represented 


by  a Fourier  sine  series 


YkC  _ Mzl 


'R  n=l 


I A.  „ sin  n 


(6.20) 
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where  A.  „ are  unknown  parameters. 

K f n 

Substitute  eq.  (6.20)  Into  eq.  (6.8)  to  get  the  final  form  of  the  down- 
wash  as 


- ? ? sin  e,. 

- - - cos  ♦^-cos 


A-l  k=1  n=1 


(6.21) 


The  tangent  flow  boundary  condition  for  thin  wings  requires  that 


(6.22) 


where  z^(x,y)  Is  the  shape  of  the  mean  camber  line.  The  N(M-1)  values  of 
A|^  ^ are  calculated  by  applying  eq.  (6.21)  at  the  chordwise  and  spanwise 
control  points  given  by  eqs.  (6.16)  and  (6.17)  and  solving  the  resulting 
matrix  equation.  For  rectangular  wings,  only  the  odd  values  of  n are 
needed  because  symmetry  properties  make  the  even  coefficients  zero. 

However,  It  Is  necessary  to  calculate  the  Aj^  ^ for  both  odd  and  even  values 
of  n for  the  skewed  wing  since  the  spanwise  variation  of  circulation  Is 
not  symmetric.  After  the  Aj^  ^ are  calculated,  the  circulation  Y|^  are 
determined  by  eq.  (6.20).  Regardless  of  the  number  (N)  of  chordwise  vor- 
tices used,  there  Is  always  a control  point  at  the  trailing  edge  which 
satisfies  the  Kutta  condition. 

The  sectional  and  wing  aerodynamic  characteristics  may  now  be  calcu- 
lated by  using  the  midpoint  trapezoidal  rule  to  reduce  the  integrals  to 
finite  sums,  as  shown  below. 


(6.23) 
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C |i^ 

*'"’j  ' 7 i *’  '**'  k-l  *’  ' “*  '''‘’ 


T 2 

St  * j ‘^m  ' *''5'  “ A J,  <'in) . ('  - »k>  »k 

-b/2  J“'  J 


(Kac/cjj  • -(yc,)^ 
Xac/c  “ "Chi/Cl 


(6.24) 

(6.25) 

(6.26) 

(6.27) 

(6.28) 


7.  LIRING  SURFACE  THEORY  FOR  SWEPT  WINGS 
The  downwash  equation  for  tapered  wings,  eq.  (6.4) , Is 


where  x = X|^^(y)  + nc(y)  . 


X-Xi 


1 + ' — - 

✓fx-x^ )^+(y-y^ )^ 


dn^  dy^  (7.1) 


(7.2) 


It  Is  Important  to  note  that  eq.  (7.1)  Is  restricted  to  planforms  with  con- 
tinuous vortex  lines.  The  geometry  of  the  swept  wing  does  not  allow  con- 
tinuous vortex  lines  In  the  spanwise  direction  because  of  a kink  at  the 
midspan,  as  shown  In  figure  3.  Therefore,  for  the  swept  wing,  eq.  (7.1) 

Is  applied  separately  to  the  right  wing  panel  and  left  wing  panel. 

Eq.  (7.1)  Is  of  the  Mangier  Integral  form  and  must  be  Integrated  by  parts 
to  obtain  the  Cauchy- type  Integral  forms  which  Is  necessary  to  apply  Lan's 
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numerical  method.  The  total  downwash  at  (n,y)  on  the  swept  wing  Is,  there* 
fore,  composed  of  the  contribution  from  the  left  and  right  wing  panels. 


w(n,y) 


!T 


Y(ni,y^) 

(y-yi ) 


7 '(yi ) 


4ir 


t,  -iU  (y-yi) 


/(x-x,)^+(y-7,)^ 


dy^  dn^ 


dy^  dn^  (7.3) 


where  x = x(n,y)  and  x^  = x.|(ni,yi).  Note  that  y^  and  y^  are  the  spanwise 
Integration  points  for  the  right  and  left  wing  panels,  respectively,  and 
y Is  the  spanwise  position  of  a control  point.  Due  to  symmetry,  control 
points  are  needed  on  one  wing  panel  only,  and  the  right  wing  panel  Is  used 
here. 

Considering  the  Integrals  In  eq.  (7.3)  separately,  the  first  Integral 
Is  Integrated  by  parts  In  the  y^ -direction,  holding  n constant,  with  y^  > 0 
and  y > 0,  as  follows: 


where 


b/2  b/2  b/2 


u E Yc  and  du 


3(yc) 

ayi 

^1 


dyi 


dv  = 


(y-yi )‘ 


1 + 


x-x. 


/(x-x^)^+(y-yi)‘ 


dy. 


(7.4) 


(7.5) 


(7.6) 


The  Intermediate  parameters  are  defined  as 


Y » d(y-y^)^  + e(y-y^)  + f 

(7.8) 

Yq  = dy2  + ey  + f 

(7.9) 

d = + 1 

(7.10) 

e - 2AB 

(7.11) 

f = 

(7.12) 

A “ ^ Alt) 

(7.13) 

B = (n-ni)(Cp  + yAm)) 

(7.14) 

"’lE  * ^LE 

(7.15) 

Am  * tan  - tan 

(7.16) 

Therefore,  the  spanwise  part  of  the  first  integral  is 


1 y(  n-j  ,yi ) c(yi ) 

4ir  y^Tj 


1 


+ 


B 


b/2 

0 


3(yc) 

3yi 


_L 

y-yi 


(7.17) 


The  circulation  y at  the  wing  tip,  y^  = b/2,  is  zero.  Integrate  eq.  (7.17) 
in  the  chordwise  direction  to  obtain 
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w^c(^*y)  + Wp(ii,y) 


(7.18) 


where  W|^^  and  W|^  are  the  downwash  components  for  the  right  wing  panel 
corresponding  to  the  first  and  second  Integrals  on  the  right  side  of 
eq.  (7.18).  The  term  Wj^  Is  the  downwash  from  the  vortices  distributed 
over  the  right  wing  panel,  while  the  term  Wj^j.  Is  the  downwash  from  a 
vortex  filament  positioned  at  mid-span  (y  = 0). 

Now,  Integrate  the  second  Integral  term  In  eq.  (7.3)  by  parts  In 
the  y^-d1rect1on,  holding  n constant,  with  y ^ < 0 and  y > 0,  as  follows 


0 

I udv 
-b/2 


0 

I vdu 

-bjz 


(7.19) 


(7.20) 


The  Intermediate  parameters  are  defined  as 

Y'  =•  d'(y-y^)^  + e’(y-7)  + f (7.21) 

Yq  » d'y^  + e'y  + f 


(7.22) 
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(A')2  + 1 

(7.23) 

2A'B' 

(7.24) 

(B')2 

(7.25) 

“ * ^LE 

(7.26) 

[2m|^g  + Am(n+ni)]  y + (n-n-i) 

(7.27) 

Therefore,  the  spanwise  part  of  the  second  integral  term  is 


0 

/ “‘''"TT 

-b/2 


Y(ni,y^)  c(y^) 


(7.28) 


Again,  the  circulation  at  the  wing  tip,  y^  = -b/2,  is  y = 0.  Integrate 
eq.  (7.28)  in  the  chordwise  direction  to  obtain 


WLc('i.y)  + W|_(n,y) 


Y(n^,0) 

y 


c 


R 


0 

i 


ayi 


y-7- 


1 L 


1 + 


dy^  dn^ 


(7.29) 


where  Wj^^,  and  W|^  are  the  downwash  components  for  the  left  wing  panel  corres- 
ponding to  the  first  and  second  integrals  on  the  right  side  of  eq.  (7.29). 
The  term  W|^  is  the  downwash  from  the  vortices  distributed  over  the  left 
wing  panel,  while  the  term  the  downwash  from  a vortex  filament 

positioned  at  mid  span. 


24 


The  total  downwash  at  (n.y)  for  the  swept  wing  Is  now  determined 
by  substituting  eqs.  (7.18)  and  (7.29)  Into  eq.  (7.3): 

w(n.y)  ■ Wpj.(n,y)  + Wp(n,y)  + WL(s(n,y)  + W|^(n,y)  (7.30) 

It  should  be  noted  that  the  two  Integrals  at  mid  span  have  the  property 
that 

(Wrc  + W|_c)  - - . 

This  Infinite  downwash  at  the  mid-span  of  a swept  wing  Is  a result  of 
the  kink  In  the  vortex  lines  across  the  center  line.  Treatment  of  this 
problem  as  pertaining  to  the  present  method  will  be  discussed  In  Chapter  9. 
The  solution  for  the  swept  wing  was  performed  by  two  different  methods 
which  are  described  In  Chapters  8 and  9. 

8.  SWEPT  WING  SOLUTION  WITH  SERIES  FOR  y DISTRIBUTION 

Similar  to  the  skewed  wing,  the  spanwise  circulation  distribution  is 
represented  by  a Fourier  sine  series.  The  circulation  at  the  chordwise 
position  X|^  = I*  ®k^  defined  as 


where y-j  * - cos  (<•*,  and  from  Appendix  A,  for  the  right  wing  panel, 

= cos"^  [|  (cos  (|)^  - 1)]  (8.2) 

R 

and  for  the  left  wing  panel 

= cos"^  [J  (1  - cos  4^) 


(8.3) 
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1. 


^ 1. 


and  are  the  unknown  coefficients. 


Defining  h tan  A^g.  m^g  s tan  A^g,  Am  = m^g  - m^g  » . 


c«  = 


_^T 


R - ^ I and  X 5 ^ , consider  each  wing  panel  separately.  The 

R 

downwash  Integrals  for  the  right  wing  panel  are  given  by  eg.  (7.18). 
Substitute  eqs.  (8.1),  (8.2)  and  the  transformed  variables  (see  Appendix  A) 


- J (1  - cos  ♦^) 

y^  > 0 

(8.4) 

* 1 (1  - cos  ♦) 

y > 0 

(8.5) 

= ^ (1  - cos  e^) 

(8.6) 

= ^ (1  - cos  e) 

(8.7) 

Into  eq.  (7.18)  and  apply  the  midpoint  trapezoidal  rule  to  get  the  downwash 
from  the  right  wing  panel  as 


N M-1  A, 


“rc^  ,j  ARTT+iTTr 


\n  [t] 

sin  e,, 

(1-cos  «j 

L.) 

'RC 


Ijk 


(8.8) 
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with  K. 


n 2 z — Z 

/t  (A,+1)(cos  (1.0-cos  (J..)  +A,B,(cos  (j.,-COS  (j.J+B, 

5 1 + — i - - ' ' i ^ (8.11) 


A^  • mj_E  + ^ (1  - cos  0,^) 

Bl  “ ®k  ■ ®i^  [m(T+aT  * ^ 


+ ^ (1  - cos 


(8.12) 


(8.13) 


The  downwash  Integrals  for  the  left  wing  panel  are  given  In  eq.  (7.29). 
Substitute  eqs.  (8.1),  (8.3),  (8.5)-(8.7),  and  the  transformed  variable 
(see  Appendix  A) 


^ (cos  (^^  - 1)  < 0 


(8.14) 


Into  eq.  (7.29)  and  apply  the  midpoint  trapezoidal  rule  to  get  the  downwash 
from  the  left  wing  panel 


"‘■Ci.j  ' wiw  J,  J, 


N M-1  A|,_-s1n  ^ .sin  0. 
r r kn  I 2 1 k 


■cos  (^, 


(8.15) 


where  K, 


/I  (A^+l)(  1-cos  ^jO^+AgBgd-cos  .frj)+B| 

= 1 ^ r» 


(8.16) 


M N M-1  s1n  Oj^-sin 

" 2AR(l+X)r(M  (2-cos  4>^-cos  ^.)  sin  4,*^ 


K,  (8.17) 

H jkt 


(Ag+l)  (2-cos  <|i^-cos  (Jij)  +A2B2( 2-cos  4)  -cos  (jijJ+B 


Ag  * - [m|_2  + Y (1-cos  0,^)] 


Am  2 

Bg  = ^ (2-cos  0,^-cos  0^)](l-cos  4)^)+  (cos  0|^-cos  0^.) 


The  control  points  are  located  at 


0^  = Y ’ ^ ^ 


= ^ . j = 1,  ....  M-1 


(chordwise) 


(spanwise) 


and  the  "joint"  In  the  horseshoe  vortices  are  located  at 


®k  “ ^~2'n  ^ ^ ^ (chordwise) 


, A = 1 , . . . , M (spanwise) 


The  downwash  equation  for  the  entire  swept  wing  is  then 
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Using  the  tangent  flow  boundary  condition  for  thin  wings  given  In  eq.  (6.22), 
N(M-l)  values  of  are  calculated  by  solving  the  matrix  equation  formed  by 
applying  eq.  (6.22)  for  the  chordwise  and  spanwise  control  points  given  by 
eqs.  (8.21)  and  (8.22).  It  Is  worth  noting  that  there  Is  always  a control 
point  at  the  trailing  edge  which  satisfies  the  Kutta  condition.  The  sec- 
tional and  wing  aerodynamic  coefficients  can  be  calculated  by  using  the 
midpoint  trapezoidal  rule  to  reduce  the  Integrals  to  finite  sums  as  shown 
In  Chapter  9. 


9.  SWEPT  WING  SOLUTION  BY  DIRECT  MATRIX  SOLUTION  FOR  Y 


As  will  be  discussed  In  the  results,  the  Fourier  sine  series  for  yc/Cj^ 
has  Instabilities  In  the  solution  which  result  from  an  Infinite  downwash 
at  the  wing  centerline.  To  obtain  a more  stable  solution,  an  alternate 
method  Is  developed  which  solves  directly  for  the  spanwise  derivative 
of  circulation  at  the  Integration  points.  Thus,  the  unknown  parameters 
In  the  matrix  equation  are  defined  as 


(9.1) 


and  similar  to  eq.  (5.16)  the  circulation  at  midspan  1s 

ir  ” 

^k,<^=0  * ■ S \,ji 


(9.2) 


The  matrix  equation  will  now  contain  N(M-l)  equations  with  NM  unknowns.  With 

N(M-l)  control  points  of  eqs.  (8.21)  and  (8.22),  N more  control  points  are 

needed  to  obtain  a solution.  It  Is  known  that  Is  finite  at  y = 0, 

^1 


thus  since  y^ 


J (1  - cos 
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= litn  "5*"  aTp  “ 1^’’'  4a“  * T5 ^ ” finite  (9.3) 

y=0  - 0 <^1  - 0 I Sin 


therefore,  = 0 which  gives  N more  boundary  conditions. 

^ 4.^=0 

Define  = tan  ■ tan  ^ “ ^TE  ~ ^LE  ~ A'r| A+1  j * 

b 

“ AR  ’ ^ ~ Paiifil  considered  separately.  The 

R 

Cauchy  singularity  resulting  from  a control  point  being  applied  on  the 
wing  centerline  is  treated  in  Appendix  B.  The  downwash  integrals  for 
the  right  wing  panel  are  given  in  eq.  (7.18).  Substitute  eqs.  (8.5)-(8.7) 
and  eq.  (9.2)  into  the  right  wing  center  integral  term,  and  replace  the 
integral  with  the  midpoint  trapezoidal -rule  sum  to  get 


'<=1.0 


N H A|,^  Sin 

k.l  «Ji  TT^oslJI 


C9.4) 


where  Kn 


^ ^ J\  (A^+l)(  1-cos  <t>j)+A^B^  (1-cos  <t»j)'*'B^ 


(9.5) 


and  A.| , B.|  are  defined  in  eqs.  (8.12)  and  (8.13),  respectively.  Substitute 
eqs.  (8.4)-(8.7)  and  eq.  (9.1)  into  the  right  wing  panel  integral  term, 
and  replace  the  integrals  with  the  midpoint  trapezoidal -rule  sum  to  get 


M N A|,,  sin  6,^ 
t=1  kil  ■ 'ijkk 


(9.6) 
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where  K. 


Mjki 


/I  (a|-H)(cos  (|ij^-cos  (>^)^+AiB^  (cos  ♦j^-cos 

" ^ ^ Bi 


(9.7) 


and  A^,  are  defined  in  eqs.  (8.12)  and  (8.13),  respectively. 


Considering  the  downwash  terms  of  the  left  wing,  eq.  (7.29),  substi- 
tute eqs.  (8.5)-(8.7),  and  eq.  (9.2)  into  the  left  wing  center  integral 
term,  and  replace  the  integral  with  the  midpoint  trapezoidal -rule  sum  to 
get 


N 

j " AR(i+X)NM 


M 


I 

l=^ 


(1-cos  (|)j)  ' 


(9.8) 


where  K, 


LC 


^ ^ /j-  (A2+1)(1-cos  »j)^+A2B2(l-cos  <t>j)+Bg 


ijk 


(9.9) 


and  A2,  B2  are  defined  in  eqs.  (8.19)  and  (8.20),  respectively.  Substitute 
eqs.  (8.5)-(8.7),  eq.  (8.14),  and  eq.  (9.1)  into  the  left  wing  panel  inte- 
gral term,  and  replace  the  integrals  with  the  midpoint  trapezoidal-rule 
sum  to  get 


M N 

'jRTPnNR  J, 


(2-cos  ♦^-cos  tj) 


(9.10) 


where  K, 


^ /^Ag+l) (2-cos  (ti^-cos  (j>j)^->-A2B2( 2-cos  <tij^-cos  ’I'jJ+Bg 


■i  jk£ 


(9.11) 
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and  A2«  B2  are  defined  In  eqs.  (8.19)  and  (8.20),  respectively. 


The  control  points  are  now  located  at 


^ . 1 = 1,  ....  N 


(chordwise) 


and 


^ . j = 0.  ....  M-1 


(spanwise) 

and  the  "joint"  In  the  horseshoe  vortices  are  located  at 


(9.12) 


(9.13) 


e - (2k-l)^  k _ 1 

2N  * ^ f • • • * 


N (chordwise) 


(9.14) 


and 


, - (2&-1)tt 

h 2M 


= 1,  ....  M 


(spanwise) 


(9.15) 


The  downwash  equation  for  the  entire  swept  wing  Is  then 


^S,j  ^.j  ‘-S.j  h.j 


(9.16) 


Using  the  tangent  flow  boundary  condition  for  thin  wings,  given  In  eq.  (6.22), 

NCM-1)  equations  are  generated  for  control  points  corresponding  to 

1 = 1,  ...,  N and  j * 1 M - 1.  When  j = 0,  eqs.  (9.4)  and  (9.8) 

are  singular.  The  limit  of  (Wpj.  + w^^^.  ) as  (|ij  -►  0 Is  approximated  by, 

1 »J  1 1 j 

(see  Appendix  B) 
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1,j=0 


^ *1  ®k 


^l“Cos  e^)) 


cos  Oj^-cos  0^ 


- Am 


(9.17) 
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+1  if  (cos  e^-cos  e,-)  > 0 

V K 1 

where  * ' 

. -1  if  (cos  e|(-cos  e^)  < 0 


(9.18) 


This  limit  is  actually  infinite  when  Am  0 and  ■Ye^'o.  and  zero  when 
Am  = 0 and  mj^^  = 0.  Equation  (9.17)  generates  N more  equations  if  applied 
when  i = 1,  ....  N and  j = 0.  Thus,  having  NM  equations  and  NM  unknowns, 
the  are  calculated  by  solving  the  matrix  of  NM  equations.  It  is  again 
worth  noting  that  there  is  always  a control  point  at  the  trailing  edge 
which  satisfies  the  Kutta  condition. 

The  sectional  and  wing  aerodynamic  characteristics  may  now  be  cal- 
culated by  using  the  midpoint  trapezoidal  rule  to  reduce  the  integrals  to 
finite  sums  as  shown  below  (for  swept  wings) 
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(9.23) 


10.  RESULTS  AND  DISCUSSION 

10.1  Prandtl's  Lifting  Line  Equation 

Figure  4 illustrates  the  convergence  of  this  method  compared  with 

the  Multhopp  series  (Fourier  sine  series)  method  of  reference  7 and  the 

conventional  VLM  (ref,  3)  for  C,  of  a rectangular  planforms  with  an  as- 

a 

pect  ratio  of  2ir.  Prandtl’s  lifting  line  theory  requires  trailing  vor- 
tices and  spanwise  control  points,  but  no  chordwise  control  points  are 
needed  because  the  theory  assumes  the  downwash  is  constant  in  the  chordwise 
direction.  Therefore,  the  control  points  are  placed  on  the  "bound"  vortex 
for  both  methods.  In  this  way,  the  accuracy  of  the  spanwise  vortex  arrange- 
ment can  be  observed  without  the  influence  of  the  location  of  chordwise 
control  points.  Figure  4 shows  that  the  present  method  converges  very 
quickly  and  approaches  the  correct  limit  as  does  the  Multhopp  series 
method.  Figure  4 shows  that  the  conventional  VLM  converges  slowly  and  c, 
does  not  appear  to  approach  the  correct  limit  when  the  curve  is  extrapolated 
to  an  infinite  number  of  trailing  vortices,  M -*■  » (1/M  -*•  0). 


I 


34 

Figure  5 shows  the  distribution  of  circulation  over  the  span  for  a 
variation  of  the  number  of  trailing  vortices.  The  circulation  distribution 
remains  virtually  unchanged  as  the  number  of  vortices  increases  from  M = 2 
to  M = 10. 

1 

10.2  Skewed  Wing  I 

t 

Table  I gives  a comparison  of  the  results  of  the  present  method  with  | 

the  conventional  VLM,  (ref.  3),  for  two  skewed  wings  with  = 30®,  I 

X > 1.0  and  AR  » 1.0  and  1.5.  The  vortex  lattice  method  is  applied  for 
4 chordwise  and  12  spanwise  integration  points.  The  results  for  the  VLM 
are  linearly  interpolated  to  correspond  with  the  present  method  results. 

This  table  shows  that  the  spanwise  variation  of  the  sectional  life  co- 
efficient compares  reasonably  well  with  the  VLM.  The  spanwise  lift  ; 

distribution  is  illustrated  in  Figure  6 for  AR  = 1.5.  The  distribution 
is  non-symmetrical  which  is  characteristic  of  skewed  wings. 

The  lift-curve  slope  is  slightly  lower  than  the  VLM  although  the  two 
are  in  good  agreement.  The  moment  coefficient  and  aerodynamic  center  j 

location  compares  reasonably  well  with  the  VLM  for  the  AR  = 1.5  wing, 
but  the  comparison  is  poor  for  the  AR  = 1.0  wing. 

10.3  Swept  Wing 

The  downwash  equation  for  the  swept  wing  yields  an  infinitely  large 
downwash  at  any  spanwise  location  where  a discontinuity,  or  kink,  exists 
in  the  spanwise  (bound)  vortex  lines  (ref.  11).  For  the  linear  taper 
swept  wing,  a kink  exists  at  the  centerline  ( y * 0)  resulting  in  an  in- 
finitely large  dowewash  at  y * 0,  which  is  unrealistic.  This  infinite 
downwash  at  the  centerline  causes  computational  problems  with  the  direct-  1 
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matrix  method  because  control  points  are  located  at  y » 0.  As  discussed 
In  ref.  11,  the  mathematical  model  of  eq.  (9.1)  Is  not  accurate  at  the 
midspan  of  the  swept  wing.  One  way  to  resolve  this  problem  too  Is  to  round 
the  leading  and  trailing  edges  near  the  midspan  point.  However,  the 
Infinite  downwash  does  not  appear  directly  with  the  Fourier  sine  series 
method  since  no  control  points  lie  on  the  wing  centerline.  On  the  other 
hand.  It  does  cause  serious  Instabilities  In  the  solution  as  the  number 
of  spanwise  control  points  Increase  because  control  points  are  approaching 
y » 0.  The  rectangular  wing  presents  no  problems  for  the  method  with 
respect  to  the  Infinite  downwash  at  y = 0 since  no  kinks  exist  In  the 
spanwise  vortex-lines  and  hence  no  Infinite  downwash. 

Table  2(a)  and  Figure  7 give  a comparison  of  the  results  of  the 
Fourier  sine-series  method  to  those  of  the  direct-matrix  method  for  a 
45®  swept  wing  of  AR  = 2.0  and  A = 1.0.  Figure  7 shows  an  Increasing 
oscillation  of  the  sectional  lift  coefficient  In  the  spanwise  direction 
for  the  Fourier  sine-series  method.  This  oscillation  Is  attributed  to 
the  control  points  near  the  centerline  of  the  wing  being  affected  by  the 
infinitely  large  downwash  at  the  centerline,  thus  Introducing  a numerical 
Instability  into  the  solution.  It  is  significant  to  note  that  the  direct- 
matrix  method  yields  a very  reasonable  spanwise  variation  of  sectional 
lift  coefficient.  Table  2(a)  shows  that  the  lift  curve  slope  for  the 
sine-series  method  and  the  direct-matrix  method  compare  very  closely, 
but  both  values  are  slightly  lower  than  the  lift  curve  slope  calculated 
by  the  VLM  (ref.  3). 

Table  2(b)  gives  a comparison  between  the  lift  curve  slopes  for  the 
direct-matrix  method  and  the  VLM  for  a 45®  swept  wing  with  taper  ratios 
of  0.0  and  0.5  and  aspect  ratios  of  2.0,  4.5,  and  7.0.  The  present  method 
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gives  very  unsatisfactory  results  for  for  swept  wings  with  taper. 

The  error  In  Increases  as  the  taper  ratio  decreased  from  1 to  0. 
a 

The  error  Is  attributed  to  the  Infinite  downwash  at  midspan,  and  a 
satisfactory  solution  would  require  rounding  the  leading  and  trailing 
edges  near  midspan  as  done  In  reference  11. 


10.4  Rectangular  Wing 

The  terms  In  the  downwash  equation  which  causes  the  Infinite  down- 
wash  at  y « 0 for  the  swept  wing  vanish  In  the  case  of  constant  taper 
and  zero  sweep.  This  result  Is  expected  since  the  rectangular  wing  has 
no  kinks  In  the  vortex  lines.  Table  3 presents  the  results  for  a 
rectangular  wing,  AR  = 2.  Table  3(a)  gives  a comparison  of  the  results 
of  the  Fourier  sine-series  method  to  those  of  the  direct-matrix  method 
for  the  variation  In  the  spanwise  sectional  lift  coefficient.  The  two 
methods  are  In  very  close  agreement.  This  adds  basis  to  the  hypothesis 
that  the  spanwise  oscillation  of  the  swept  wing  sectional  11ft  coefficient 
using  the  Fourier  sine-series  method  Is  attributed  to  the  infinite  down- 
wash  at  y = 0,  which  causes  numerical  Instability  In  the  solution. 

Table  3(b)  shows  that  the  direct-matrix  method  yields  a variation  of 
sectional  lift  coefficient  very  close  to  that  of  reference  7.  Some  dif- 
ference In  the  two  methods  may  be  as  a result  of  the  quadratic  Interpola- 
tion method  used  to  obtain  the  comparable  sectional  lift  coefficients. 
Table  3 also  shows  that  the  direct  matrix  method  gives  the  same  overall 
aerodynamic  coefficients  as  those  In  reference  7.  The  results  for  the 
overall  aerodynamic  coefficients  for  the  Fourier  sine-series  method  do  not 
compare  as  well  as  the  direct  matrix  method.  The  major  difference  between 
the  present  methods  and  the  method  of  ref.  7 Is  that  the  present  methods 
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apply  the  semi -circle  distribution  of  control  points  and  Integration 
points  over  the  two  half  spans  and  the  method  of  ref.  7 applies  the  semi- 
circle distribution  over  the  whole  span. 

10.5  Compressibility  Effects 

All  of  the  results  given  In  this  report  are  for  Incompressible  flow. 
Subsonic  compressibility  effects  may  be  Included  very  simply  by  applying 
the  Prandtl-Glauret  rule  (ref.  3). 

10.6  Computational  Requirements 

The  examples  computed  herein  were  all  performed  on  an  IBM  370/165 
computer.  Typical  cases  for  the  lifting  surface  theory  required  about 
20  seconds  of  CPU  time. 


11.  CONCLUSIONS 

The  following  conclusions  are  drawn  from  the  present  Investigation; 

1.  The  new  lifting-surface  method  developed  herein  was  found  to  compare 
well  with  other  lifting-surface  theories,  but  with  much  smaller 
computational  times. 

2.  This  new  method  Is  more  accurate  and  converges  faster  than  conven- 
tional vortex  lattice  methods. 

3.  For  skewed  and  rectangular  wrings,  the  results  using  the  Fourier-sine 
series  for  the  spanwise  variation  of  circulation  are  essentially  the 
same  as  those  using  the  direct  matrix  solution. 

4.  For  swept  wings,  considerable  differences  were  found  in  the  results 
using  the  Fourier-sine  series  as  compared  to  the  direct  matrix  solution. 
The  Fourier-sine  series  method  was  found  to  give  an  oscillatory  solution 
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for  the  circulation,  and  the  amplitude  Increased  as  more  terms  were 
Included  In  the  series.  The  direct  matrix  solution  was  found  to 
give  more  satisfactory  results. 

The  basic  equation  used  In  lifting  surface  theories  gives  an  Infinite 
downwash  at  the  midspan  of  swept  wings.  This  singularity  In  the 

equation  Is  caused  by  the  "kink"  at  midspan  In  the  spanwise  vortices, 
and  the  effects  of  It  can  be  minimized  by  positioning  control  points 
away  from  the  midspan  position.  Another  alternative  Is  to  round 
the  leading  and  trailing  edges  near  midspan,  thereby  eliminating  the 
kinks  In  the  spanwise  vortices. 
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APPENDIX  A 

Transformations  for  Swept  Wings 

The  present  method  uses  the  "semi-circle"  coordinate  system  for  the 
spanwise  and  chordwise  direction.  The  "semi-circle"  method  Is  applied 
separately  to  the  right  and  left  wing  panels  for  the  swept  wing  because 
of  a kink  at  the  midspan,  (see  figure  A-1).  The  transformation  from 
(n,y)  to  (0,(^)  Is  made  by  Introducing  the  transformations: 

Chordwise  Direction 

n ® ^1  - cos  e)  control  point 

^1  - cos  e^)  vortex  point 

where 

e » 0 at  T)  = 0 (leading  edge) 

0 = ir  at  n = 1 (trailing  edge) 

Spanwise  Direction 

Right-Wing  Panel;  y = ^(1  - cos  4.)  control  point  (A-3) 

y-j  ' ^1  - cos  4.^)  vortex  point  (A-4) 

where 

(>  = 0 at  y ■ 0 

♦ » IT  at  y “ ^ 


(A-1) 

(A-2) 


J 
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Left-Wing  Panel:  = ^cos  4i  - 1)  vortex  point  (A-5) 

where 

■ 0 at  y ^ • 0 

♦ ^ = IT  at  y-|  = - I 

The  control  points  are  only  applied  to  the  right  wing  panel  because  of 
syimietry,  thus,  eq.  (A-3)  applies  for  the  left  wing  panel  integration 
also. 

For  the  case  where  circulation  is  represented  by  a Fourier  sine 
series,  it  is  necessary  to  transform  the  half-span  "semi-circle"  coordinate 
system  into  a full  span  "semi-circle"  coordinate  system.  This  is  accom- 
plished by  defining  for  the  right  wing  panel 


^ 11 

= cos"  [^cos  - 1)]  (A-6) 

R 

where 


(j*!^  = ^ at  = 0,  y = 0 


and  for  the  left  wing  panel 


where 


4*^  “ cos"^  [^1  - cos 


= J at  » 0,  y = 0 

■ 0 at  ?i  “ If,  y * - |- 


(A-7) 


APPENDIX  B 

Singularities  at  Midspan  of  Swept  Wings 

Singularities  occur  In  eqs.  (6.21)  and  (6.33)  when  y y -►y^. 
and  y -»■  0.  These  Cauchy  singularities  are  treated  using  Lan's  approach 
(ref.  4).  The  downwash  equation  for  the  right  wing  panel  In  Integra! 
form  Is 


VrRsil  s1n9,-l 


cos  di,-COS  ^ 


COS  ill. -COS  (ji 


COS  4. -cos  4 


I 


h 


, «,  = 1 , . . . , M 


(spanwise) 


M ’ 


j = 0,1,  ...,  M - 1 (spanwise) 


It  Is  known  that  is  finite  at  y = 0,  thus  since  ^ (1-cos 


(l^l 


3v 

11m 


1^  h c?n  ^ 
*^-.0  ^*1  ^ *1 


therefore 


This  shows  the  extreme  right  side  term  of  eq.  (B-5)  to  be  zero  when  j = 
thus  eq.  (B-5)  becomes 


ICfjJ  ^n('fr*^l)d^i  _ M 3(^,  R^^Lo 


COS  (^,-COS  t 


cos  ♦»-cos  * 


(A,+l)(cos  ♦ -cos  +A,B,(cos  *.-cos  *4)+B 


Therefore,  replacing  eq.  (B-1)  by  the  midpoint  trapezoidal -rule  sum  to  get 
the  downwash  from  the  right  wing  panel 
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Ji=l  k=l 


y 

cos  *^-cos 


(B-11) 


where 


(B-12) 


The  significant  result  is  that  eq.  (B-11)  applies  for  control  points 

j » 0 M - 1,  giving  NM  equations  with  NM  unknowns. 

It  can  be  shown  similarly  by  using  Lan's  method  for  treating  Cauchy 
singularities  that  the  downwash  from  the  left  wing  panel  is 


M N Sin  0^ 


„ _ IT  Y Y KX. K ^ 

\ ^ j ■ AR(1+X)NM  (2-cos  .^^-cos  ' \ 


ijkt 


(B-13) 


where 


‘ijkt 


y^(A|+l) (2-cos  4'j^-cos  (j).)^+A2B2( 2-cos  (|>j^-cos  i-J+B 
= T + J 


(B-14) 


and  the  integration  points  and  control  points  are  given  in  eqs.  (B-6)  and 
(B-7)  respectively. 

The  center  integral  terms  of  eqs.  (6.21)  and  (6.33)  have  a singularity 
as  y -*•  0.  Add  the  center  integral  terms  to  get 


I Y(n^.O) 

y 

0 

“T‘ 

“lc  * "rc  ■ ^ 


(B-15) 


1 


T 
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The  limit  of  this  equation  as  y 0 is 


1 


h-ni| 


in-ni  I In  n-i  I 

T(n,.0)  7 ^ 

I in-Oi/ 


dn,  (B-16) 


where  n = ^1-cos  e),  dn^  = j sin  de^ 


Consider  the  integral 


r ln-ni | r |cos  e,-cos  0| 

' y dn-,  = 5-  sin  e,  de, 

^ (cos  e^-cos  e)^  ' ' 


0 


lim 


e 0 


e-e 


sin  9 


1 


cos  9.J-C0S  0 


d0 
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I 


sin  0^  d0^ 


0+e 


COS  0-|-cos  0 


lim 
e ->  0 


|Tog  sin  (0  - je)  - log  sin  ]-e 


" sTnTT 


7 


(B-17) 


1 


!l 


II 


Thus,  the  downwash  at  the  center  y » 0 is  infinite  for  a swept  or  tapered 
wing.  Note,  however,  that  when  lUj^^  = 0 and  Am  = 0,  i.e.,  a rectangular 
wing,  the  limit  of  eq.  (B-16)  is  zero. 


For  computational  purposes,  eq.  (B-16)  is  replaced  by  the  midpoint 


trapezoidal  rule  sum  to  get 


j=0 


M N 

I I A, 

k=l 


U ^1  ®k 


2(mL£+  ^1-cos  0^)) 


cos  Sj^-COS  0^ 


- Am 


Mhere 


r+1 


>1 


If  (cos  0J,  - cos  0^)  > 0 
If  (cos  0,^  - cos  0.. ) < 0 


(B 


(B 
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TABLE  1.  Results  for  Skewed  Planforms 
X=1.0,  AR»1.5  and  1.0 


Values  of  C /C, 
A L 


.9659 
.8660 
.7071 
0.5000 
0.2588 
0 

-0.2588 

-0.5000 

-0.7071 

-0.8660 

-0.9659 


Present 
N*4,  M=12 


1.3517 
1.6742 
1.9416 
.1339 
.2408 
1.2620 
1.2032 
1.0718 
0.8748 
.6209 
.3230 


VLM 

(ref.  3) 
N=6,  M=40 


1.4032 
1.6908 
0.9423 
1.1258 
1.2276 
1.2464 
1.2003 
.0608 
1.8719 
i.6331 
.3682 


.9659 

.8660 

.7071 

.5000 

.2588 

0 

.2588 

.5000 

.7071 

.8660 

.9659 


Present 
N=4,  M=12 

VLM 

(ref.  3) 
N=6,  M*40 

m 

0.3884 

0.6683 

0.9184 

0.9185 

1.1172 

1.1085 

1.2363 

1.2227 

1.2698 

1.2544 

1.2186 

1.2042 

1.0880 

1.0778 

0.8850 

0.6406 

0.3252 

0.3715 

Overall  Values 

Present 

VLM 

Present 

^L 

1.9629 

2.0057 

a 

”^M 

0.4120 

0.4417 

a 

Xac/c 

0.2048 

0.2202 

^L 

1.4399 

1.4736 

a 

0.2581 

0.1546 

a 

Xac/c 

0.1793 

0.1049 

TABLE  2.  Results  for  45®  Swept  Wing 
AR«2.0,  X^l.O 

Values  of  C^/Cj^ 

N=4,  M=12 


Sine  Series  Non-Series 
Method  Method 


.0170 
.0670 
.1464 
.2500 
.3706 
.5000 
.6294 
.7500 
.8536 
.9330 
0.9830 


.1515 

.1595 

.1699 

.1794 

1.1547 

1.1557 

.9611 

.0740 

0.4839 

0.7937 

0.0749 


.1550 
.1623 
.1719 
1.1775 
.1660 
.1250 
.0430 
.9130 
.7349 
0.5152 
0.2655 


Overall  Values  of  C 


Sine  Series 
Method 

Direct  Matrix 
Method 

VLM 

(ref.  3) 

2.2544 

2.2586 

2.265 

Swept,  Tapered  Wing 
A. .=45° 


N»4  M=12 


Direct  Matrix 
Method 

VLM 

(ref.  3) 

2.6463 

2.4408 

3.6774 

3.3690 

4.0393 

3.7013 

9.3418 

2.3606 

11.7134 

3.4779 

f 


I 

I 
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TABLE  3.  Results  for  Rectangular  Wing 
AR-2.0 


Values  of 


a)  Comparison  of  Fourier  Sine-Series 
Method  to  Direct  Matrix  Method 

N-4,  M»8 


2y/b 

Sine  Series 
Method 

Direct  Matrix 
Method 

0.0381 

MEM 

MSM 

0.1464 

mEBM 

0.3087 

1.2023 

0.5000 

1.1059 

0.6913 

0.9324 

0.9358 

0.8536 

0.6820 

0.6850 

0.9619 

0.3618 

0.3636 

b)  Comparison  of  Direct  Matrix  Method 
to  Method  of  Reference  7 

N»4,  M=8 


2y/b 

Direct  Matrix 
Method 

(Interpolated) 

Method  of 
ref.  7 

■H 

mSM 

IHH 

mm 

msM 

IM 

1.0651 

0.7071 

0.9171 

0.8315 

0.7237 

0.7290 

0.9239 

0.4997 

0.5084 

0.9808 

0.2575 

0.2618 

Overall  Values 


Sine  Series 
Method 

D1 rect-Matrl x 
Method 

Method  of 
ref.  7 

2.4733 

2.4732 

2.4732 

a 

o 

0.5143 

0.5187 

0.5187 

0.2079 

0.2097 

0.2097 

x/c  ■ (1  - cos  0)/2 


0 Vortex  Position:  (k  » 1,  ..., 

X Control  Points  ; (i  = 1,  ...» 

Figure  1.  Lan's  vortex  arrangement  for  airfoils 


o VLM  (ref.  3) 

□ Fourier-Sine-Series  (ref.  7) 
A Present  Method 
M No.  of  Spanwise  Control  Pts. 


.8  L/ 


FIGURE  4.  Convergence  of  C^^  Predicted  by  Three  Methods  of  Solving 
Prandtl's  Lifting-Line  Equation,  Rectangular  Wing,  AR=2ir. 


